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Abstract 

D— dimensional central and complex potentials of a Coulomb plus quartic-polynomial form are 
considered in a VT— symmetrized radial Schrodinger equation. Arbitrarily large finite multiplcts 
of bound states are shown obtainable in an elementary form. Relations between their energies and 
couplings are determined by a finite-dimensional secular equation. The Bender's and Boettcher's 
one-dimensional quasi-exact oscillators re-emerge here as the simplest chargeless solutions. 
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1 Introduction 



Recently, Bender and Boettcher [1] studied the three-parametric family of Hamiltonians 

H {BB) = -d%- x* + 2iax 3 + cx 2 + i(a 3 - ac-2J)x, a,ceM, J = 1,2,... (1) 

defined on a certain complex curve x — x(t) S(F, t € (—00,00). Their main result was an explicit 
construction of a J— plot of bound states in an elementary and closed form 

^(x)=e-* x3 / 3 - a * 2 / 2 -* bx P jJ _ 1 (x), b={a 2 -c)/2 1 j = 1,2,..., J. 

Here, denotes certain polynomials of the (J — l)-st degree. The normalizability of the 

multiplet is manifestly guaranteed via pre-selection of the asymptotics of x(t), 

x(t) ~ (cos ip ± i sin ip) t, t — > =Foo , ip e (0, tt/3). 

It is easy to check that we may even stay on the straight line x(t) — t — is, provided only that we 
choose it in such a way that Re(2e + a) > 0. 

An older paper by Buslaev and Grecchi [2] may be recalled for the rigorous proof that the 
Hamiltonian (1) has the real and discrete spectrum. The paper shows that after a constant shift 
S ~ b 2 + J a/2 the spectrum will coincide with the energy levels of the D— dimensional anharmonic 
oscillator H {AH °^ = -A + jr 2 + \[r 2 ] 2 in its £-th partial-wave mode such that 2 J = D + 
2i — 2. This assignment (or "BG transformation") combines a change of variables with Fourier 
transformation. Its application to H^ BB ^ at an integer J specifies either our choice of the AHO 
partial wave £ = 0, 1, ... or of the (even) dimension D. In the opposite direction it throws new light 
on the non-Hermitian Hamiltonians which exhibit the puzzling [3] VT— symmetry 

VTH (BB) = H {BB) VT ^ = ^,(_ a .) j T*P(x) = i>*{ X ). 

The J-plets of their exact states become transformed in the unusual anharmonic oscillator states 
obtainable in terms of certain elementary Fourier- type integrals (cf. also rcf. [4] in this context). 
This adds a new and strong reason why the models of the type (1) deserve a thorough attention. 

We may only feel dissatisfied by a certain incompleteness of the whole picture: Why the BG 
transformation prefers the even dimensions? What could be done at the odd D = 1,3, . ..? The 
free variability of the integer D would be highly welcome also in some phenomenological D ^> 1 
models in nuclear physics [5] , quantum chemistry [6] and atomic physics [7] . 

The second reason why the class of models (1) looks so inspiring is related to the second version 
of the BG transformation, presented in paper [2] as "the second main result" . It starts from a less 
usual, volcano-shaped model V(r) ~ uP'r 2 — [r 2 ] 2 . The current partial-wave decomposition of its 
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D— dimensional wave functions 

oo 

= r^~ D ^ 2 ip{r) x angular part , r = |f| € (0, oo) 

reduces its Schrodinger equation to the ordinary equation on the half-axis. Using a complex shift 
of the coordinates again, Buslaev and Grecchi arrive at a new, properly VT~ symmetrized two- 
parametric operator 

H ( BG ) = _ d 2 + HL+1) + u 2 x 2 _ x a ; x = x (t) = t -ie, L = £+(D-3)/2, (2) 

extended to the whole real line t E M and acting in the current Hilbert space L 2 (1R). They prove 
its isospectrality with a certain one-dimensional (and Hermitian) double-well oscillator. In the 
resulting quadruple scheme 



1 — dimensional 
Vi — symmetric " 
-x 4 BB model (1) 

I BG transformation 



1 — dimensional 
and Hermitian 
+x 4 double well 



| BG transformation 



D — dimensional 
and Hermitian 
+x 4 AHO well 



x^i (t — i e) 



VT — symmetrized 

and regularized 
-x 4 "volcano" (2) 



the two columns are related by the elementary changes of variables. The upper-left-corner item is 
characterized by its partial solvability. In what follows we intend to show that the lower-right-corner- 
like singular Hamiltonians may equally well exhibit the same (or at least very similar) multiplet or 
"quasi-exact" [8] solvability. 



2 VT— regularized singular oscillators 



The problem we shall solve here is the Schrodinger equation 
dx 2 



L(L + 1) .d. ,2-i 4 
= Yi — \-icx + bx +iax ■ —x 



ip{x) = Eip(x) 



(3) 



x" x 

defined on curves x(t) ^ to be specified later. It contains the bound state problems (1) and 
(2) after an appropriate choice of its parameters. As long as certain characteristic features of the 
similar PT— symmetric systems will be of an immediate relevance here, let us first review some of 
them briefly. 
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2.1 VT— symmetry 

Many complex models with VT symmetry (and, let us assume, real spectra, Im E = [9]) are 
defined on real line [10]. Their potentials comprise a spatially symmetric real well plus its purely 
imaginary antisymmetric complement. An interpretation of these models parallels the usual real 
and symmetric bound state problems in one dimension since one may try to switch off the imaginary 
force smoothly. Typically, the spectrum splits in the convergent and divergent parts in this limit. 
Recalling cq. (2) for illustration and definiteness, we identify the convergent part with the standard 
physical solutions ip(qo)( r ) ~ r ^+( £) - 1 )/ 2 Q f a "quasi-odd" or "regular" type. The divergent part 
represents the "quasi-even" components ip( qe )(r) <~ r (3--D)/2-£_ p or a f ur ^ ner explicit illustration 
one may recall the harmonic oscillator [11] or a few other solvable examples [12]. 

The less trivial set of the models which "weaken" their Hermiticity to the mere VT symmetry 
may be defined off the real line [13]. In the usual mathematical notation one speaks about sectors 
where the L 2 (M) boundary conditions are imposed. Thus, returning to our examples (1) - (3) we 
define sectors Sk — {x e C; x ^ 0, | arg(x) — \{2k — l)ir\ < ^tt} and choose S4 and Sq in accord 
with the recommendation of rcf. [1], or Si and S3 as accepted in ref. [2] (cf. Figure 1). 

The latter ambiguity of the curve x(t) may be characterized by a "signature" a = ±1 which easily 
distinguishes between the two T— conjugate curves x(t) — x^(t) with the asymptotics bounded by 
the ip — ► 7r/3 Stokes lines, 

x^\t) ~ (cos ip ±ia sin <p)t, t — » =Foo , tp € (0, 7r/3). (4) 

Even in the admissible straight- line or tp — > extreme x^(t) —t — iae, the possible branch-point 
singularity in the origin must be avoided properly, e.g., in the trivial T— conjugation manner which 
leaves the (real) spectrum unchanged. For this reason the condition which determines e is now 
more restrictive and reads 2s > \a\. 

We may summarize that our curves of integration x{t) — x^\t) may be straight lines or 
hyperbolic shapes open downwards (<r = +1, [1]) or upwards (a = —1, [2]). For the fully regular 
forces their non-asymptotic form may be deformed almost arbitrarily. For all the other models with 
singularities these curves must avoid all the cuts. In particular, eqs. (2) and (3) admit the presence 
of an essential singularity in the origin. After we cut the plane upwards (for a = +1) or downwards 
(for a = —1), both our curves x = x^ (t) may be kept smooth and, say, symmetric with respect to 
the T reflections. 

As long as we are relaxing the usual Hermiticity ^ [H^ BB ^} + , both the two independent 

wave function solutions remain equally admissible in the vicinity of the origin. This is the main 
formal reason why the Bender's and Boettcher's complex oscillator remained solvable: One only has 
to combine the quasi-even behaviour ip( qe )(r) ~ O(r ) with its quasi-odd parallel ip( qo ){r) ~ ©(r 1 ) 
in the single ansatz. Here we intend to proceed in the same manner. 
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2.2 Recurrences 

Schrodinger equation (3) with the integer or half- integer values of L = £+ (D — 3)/2 = K/2 remains 
perfectly regular along both our non-self-intersecting integration paths. After a change of variables 
x^(t) = —iy(t) we get the new differential equation for tp(x) = x{^ x ) which contains, formally, 
no imaginary units, 

x(y) = E x (y). (5) 

We shall search for its solutions by a power-series ansatz, assuming that such a series terminates. 
Demanding also the manifest normalizability in a way depending on the signature a = ±1 wc arrive 
at a virtually unique formula 

/If \ N 

x(y)^cxp(a-y 3 + -Ty 2 + Sy) ]>>„2,"- L . (6) 

It combines both the respective y~ L and y L+1 quasi-even and quasi-odd components, and equation 
(5) determines the asymptotically correct values of our auxiliary parameters as well as the unique 
termination-compatible coupling c, 

T = a/2a, S = (b - T 2 )/2a, c = c(N) = -2TS - a ■ (2N - 2L + 2). 

The resulting ansatz (6) represents the desired normalizable solutions if and only if its coefficients 
h n are compatible with the recurrences 

h n+ iA n + h n B n + h n -iC n + h n - 2 D n = 0, n = 0, 1, . . . , N + 1. (7) 

Their coefficients 

A n = (n+l)(n-2L), B n = S(2n-2L)-d, C n = S 2 +T(2n-2L-l)-E, D n = 2a(n-N-2) 
are all elementary. 

3 Finite-dimensional Schrodinger equation 

In the light of their strict postulated termination, our recurrences (7) represent just a finite set 
of N + 2 linear algebraic equations for TV + I unknown coefficients h n . Such a set is, obviously, 
over-determined. Its (N + 2) x (N + 1)— dimensional non-square matrix has in effect two main 
diagonals, 

B n = S(2n -K)-d, C n = S 2 + T(2n - K - 1) - E, K = 2L. 



dy 2 



L(L + 1) d 

y 2 y 



cy-by -ay 3 -y 
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Both the quantities d and E can play a role of an eigenvalue simultaneously [14]. In the other words, 
we may interpret the whole linear set of equations as the two coupled square-matrix problems 
indicated, say by the single and double line in the whole non-square system 

/ B Ap 

C x B x A x 
D 2 C 2 '•• 

Bn-i An_x 
Dn Cn Bn 
\ Dn+x Cn+i j 

In both cases (i.e., omitting the last or first line, respectively) we get a routine secular determi- 
nant with an eigenvalue on its main diagonal. Both these two independent secular or termination 
conditions have to be satisfied simultaneously. 



f ho \ 

hi 
hi 

\h N ) 



= . 



(8) 



3.1 Auxiliary constraint 



As long as the upper diagonal A n = (n + l)(n — K) vanishes at n = K , there emerges an important 
asymmetry between our two eigenvalues d and E. Whenever one keeps just a few lowest partial 
waves, we may say that the integer K = 2L remains "very small" , at least in comparison with N 
which can/should be "large" or at least "arbitrary" in principle. This is our present fundamental 
observation. As a consequence of the related disappearance of Ak =0we shall achieve a thorough 
simplification of the construction of bound states. 

At the lowest possible K = we return immediately to the Bender's and Boettcher's proposal. 
Their choice of d = in cq. (1) (with J = N+ 1) gives us the highly welcome possibility of omitting 
the whole first row from eq. (8). The rest of this equation is the usual square-matrix diagonalization. 
All the exceptional quasi-exact eigenvalues Ej, j = 1, 2, . . . , N + 1 may be determined as roots of 
a polynomial of the (N + 1)— st degree [1]. 

A transition to the nonzero integers K is easy. In place of using the trivial d = 0we have to 
satisfy the (K + 1)— dimensional sub-equation 



det 



' Bo 


Ao 


Ci 


Bi 


D 2 


c 2 



\ 



Ai 



V 



D 



K 



Br-i 
C k 



A K -i 
B k ) 



0. 



(9) 



This can fix, say, the eligible electric charges d as functions of the other parameters. These functions 
can be used as a starting point of a facilitated solution of our original problem (8). 
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3.2 Closed formulae for the energies 

In the simplest test of our new and general recipe let us return, once more, to K = in eq. (9) and 
derive 

d = 0, K = 0. 

This reproduces the model H^ BB ^ of rcf. [1]. In the next (and the first really innovative) K = 1 
step we get the double root d = dr±^ (E) = ±\/E. Each choice of the sign specifies a different 
relation between the charge and the energy. Both these signs remain admitted by the inverted, 
unique recipe 

E = E K (d)=d 2 , K = l. 

It defines the energy as a function of the charge. The numbering of the separate elements of our 
new, K = 1 multiplets of bound states is in this way transferred directly to the admissible charges 
d = dj. Formally this means that the energy E may be eliminated from all our algebraic equations. 
The values of the charge d remain the only unknown quantities. 

At the subsequent integer K — 2 we get, with a bit of luck, a highly compact energy formula 
again, 

, . d 2 ST + aN 
E = E K {d) = — +2 , K = 2. 

The further step gives the two rules or roots E 3 (d) = F of the quadratic equation 

9F 2 - 10 d 2 F + d 4 +48dST + 48daN - 24 d- 72 S- 36 T 2 - 0, K = 3. 

Hence, all the K > 3 cases require the so called Grobner elimination which should (and easily can) 
be performed by a computer. We only quote a similar calculation [6] and omit the further details. 



4 Secular determinant and its roots 



In the usual quasi-exact manner [8] the multiplets of the (real) values of dj have to follow from our 
termination postulate at any positive integer TV. Let us now assume that the value of K is fixed. 
After the insertion of its energy formula E = Ex(d) in the recurrences or equation (8) we are left 
with the square-matrix secular equation of the dimension (N + 1) x (N + 1), 



/ did) Bi{d) A l 
D 2 C 2 (d) ••• 



det 



D 



N 



\ 



C N (d) B N (d) 

Dn+1 Cjy + i(d) J 



. 
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This equation has to determine the unknown eigencharges. In its analysis we may skip the already 
known Bender's and Boettcher's K = case [1] and pay attention to the next few K = 1, 2, . . .. 



4.1 The first nontrivial option: K = 1 

Insertion of the available formulae gives us the K = 1 secular determinants with the only N— and 
a— dependent elements Dj(N) = 2<r(n — N — 2), 

/ S 2 -d 2 S-d \ 

D 2 (N) S 2 + 2T-d 2 3S-d 3 

dot. D 3 (N) S 2 +AT-d 2 '■■ '■■ =0. (10) 

y D N+1 (N) S 2 + 2NT-d 2 J 

At any N this equation gives an "exceptional" explicit root d~S. Although it looks like an artifact 
of our present approach, its further inspection reveals its acceptability. For example, at N = 1 this 
root gives a solution provided that we have ST = —<r/2. 

All the other roots of eq. (10) remain manifestly N— dependent. They also exhibit certain 
symmetries. For example, the change of the signature a — ► —a may be compensated by the 
simultaneous sign-change of S — ► — S and d — > — d. 

In a more constructive setting it makes sense to pick up a specific N. Starting from the smallest 
N = 1 we arrive at the cubic equation 



-d 3 - S d 2 + (S 2 + 2T) d + 2 a + S (S 2 + 2T) = 0, 



N = 1. 



We may determine the S— and T— dependence of its three roots d\, c?2 and d% via Cardano formulae. 
If needed, the boundaries of the domain of their reality may be determined numerically, in a complete 
parallel to the K = study [1]. In order to prove that this domain is not empty at K — 1, we 
may recall a sample triplet, say, of d h2 ,3,{ S i T ) = (-5.303953910,-3.103253421,5.407207331) at 
(5,T) = (3,10). 

The growth of the dimension N + 1 makes the practical determination of the roots d of eq. 
(10) more complicated. This is well illustrated by its next two K = a = 1 explicit polynomial 
representants, 

d 5 + Sd 4 - (2S 2 + 6T) d 3 -(6 + 2S 3 + 6ST) d 2 + (6 S 2 T + 4 S + 8 T 2 + S 4 ) d 
+ 10 S 2 + 16T + 6 S 3 T + 8 ST 2 + S 5 = 0, N = 2 
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-d 7 - Sd 6 + (12 T + 3 S 2 ) d 5 + (12 + 12 ST + 3 S 3 ) d 4 - (16 S + 44 T 2 + 3 S 4 + 24 S 2 T) d 3 
- (40 S 2 + 88 T + 44 ST 2 + 3 S 5 + 24 S 3 T) d 2 
+ (12 + 16 S 3 + S 6 + 44 S 2 T 2 + 64 ST + 48 T 3 + 12 S 4 T) d 
+ 152 S 2 T + 28 S 14 + 84 S + 144 T 2 + S 7 + 44 ^T 2 + 48 ST 3 + 12 S 5 T = 0, = 3. 

To this list with the respective 16 and 31 terms one could add the next N — 4 item containing 53 
terms, etc. This would complement the similar K = formulae displayed in detail in ref. [1]. 

4.2 Numerical illustration 

An overall insight in the structure of the spectra may be mediated by their simplest S = T = 
sample. A posteriori, the main merit of such a choice may be seen in a drastic formal reduction of 
the underlying secular polynomials. One may factorize many of them by the purely non-numerical 
means. 



4.2.1 K = 

In the paper [1] the domain of existence of the full quasi-exact N— plets with K = has been 
determined numerically. This analysis excluded the point S = T = where, empirically, the VT 
symmetry becomes spontaneously broken. Computationally, this means that some of the energies 
coalesce and move in pairs off the real line. Hence, our understanding of the K = S = T = 
spectra is not complete yet. 

Direct computation at = 0, N = 1 and N = 2 gives just the trivial E = 0. In the next two 
cases we get the a— independent equations 

E 4 - 96 £ = 0, 
-E 5 + 336E 2 = 0, 

with the single nonzero real root E <~ 4.48 at A^ = 3 and E ~ 6.95 at A^ = 4. For the general N 
the secular equation reads 

/ —E 1-2 \ 
2No~ -E 2-3 



N = 3 
N = 4 



det 



V 



6cr -E (N — l)N 
4cr —E 
2d —E 



0. 



Its a— independence and polynomiality in x = E 3 ^ can be proved in an easy exercise. As a 
consequence, its roots remain non-numerical up to A^ — 13. In the latter extreme we get the real 
quadruplet of nonzero energies E - 9.381, 17.768, 26. 487 and 35.535 . 
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4.2.2 K=l 

After one moves to the negative signature a = — 1, only certain signs change in the secular polyno- 
mials of sect. 4.1. At both a = ±1 we may reduce the first few cases to their pertaining S = T = 
forms. They remain non-numerical and exactly solvable up to the dimension N = 6. Their first 
few samples are 

-d 3 + 2a = 0, A = 1 

d 5 -6crd 2 = 0, A = 2 

-d 7 + 12crd 4 + 12d=0, A = 3 

rf 9 -20crd 6 -76d 3 + 512cr = 0, N = A. 

In detail, the single real root d rts 1.26 • cr is obtained at TV = 1, and one double zero and one real 
nonzero root d w 1.71 ■ a appear at A = 2. One simple zero and one positive and one negative 
root (d\ s=s 2.35 • (7, d 2 ~ —0.975 • <r) follow at A = 3 while, finally, three nonzero real roots 
& ■ ^1,2,3 ~ (2-82, 1.55, —1.83) emerge from our last displayed polynomial at N = 4. 

4.2.3 A = 2 

In a move to the higher K > 1 one has to notice that the first few choices of A < K are too formal 
and do not make an explicit use of the present "advantage" Ak = at all. In this sense, the K = 2 
illustration has to start from the first nontrivial A = 3. After we fix a = +1 for brevity, we get 
the secular polynomial of the deterring twelfth degree. Fortunately, in the same spirit as observed 
above the abbreviation x = d 3 ^ reduces it again to a (still solvable) quartic equation, 

x 4 + 331776 - 96a; 3 + 384 x 2 + 18432 x = 0, a = 1, K = 2, A = 3. 

This equation possesses the two real and positive roots, namely, x\ — 24 and 

3 1 ^ 64 
x 2 = 24 + 16 V9 + Vfi + - w 88.87294116. 

V9 + \/T7 

These roots lead to the real and positive charges d\^ ~ (2.88, 4.46) and to the K = 2 energies 
E2(dj) as prescribed above. The parallel problem with a = — 1 leads to the full quadruplet of the 
negative real roots 

2+2,3,4 ~ (-199.78, -72.00, -14.65,-1.57), a = -1 

and to the eigencharges with the same signs, ^1,2,3,4 ~ (—5.84, —4.16, —2.45, —1.16). 

We may conclude that at a fixed K (i.e., for a specific partial wave t) we get in general a 
multiplct of states connected by some broken lines in the two-dimensional charge-energy plane. 
Only in the simplest K = special case this "guide for eye" becomes the Bender-Boettcher straight 
line d = const = 0. 
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5 Discussion 



We have shown that the N— plets of exact bound-state solutions of our general quartic problem 
(3) can be constructed in closed form at any integer degree TV and dimension D. In the other 
words, under certain relationship between the couplings and energies, arbitrarily large multiplcts 
of solutions (with real energies!) proved obtainable from a single and finite-dimensional secular 
equation. The Bender's and Boettchcr's quartic example re-appears here as the simplest one- 
dimensional special case with the unique (and, incidentally, vanishing) Coulombic charge d. 

In many applications of practical interest the closed quasi-exact solutions are of the similar 
Sturmian type. Their bound states are numbered by one of the couplings but still share, mostly, 
the same value of the energy. Many Hcrmitian models belong to this category, be it one of the most 
popular non-polynomial examples [15] or one of the historically first models of the quasi-exact type 
[16]. Only the most popular partially solvable sextic model [17] owes for its popularity to its most 
standard variable-energy character. 

We have presented the explicit construction of a few multiplets which lie along the more general 
curves (e.g., parabolas) in the energy-coupling space. In combination with their overall VT sym- 
metric quantum mechanical framework these multiplets offer several new challenges, e.g., in their 
completely missing abstract interpetation, say, from the modern Lie-algebraic point of view [18]. 

In our present approach the recent unexpected construction of Bender and Boettcher [1] finds 
one of its simplest and quite natural explanations. We have found that its solvable status has 
been mediated first of all by its complex, non-Hcrmitian background. In a way which clarifies the 
whole paradox the exceptional features of the Bender's and Boettchere's oscillator were related to 
its special chargeless form. 

In another approach to the same problem one may just speak about the underlying system 
of recurrences. In this language, it remained unnoticed in the current literature on anharmonic 
oscillators that the coefficient Ak can vanish at the integer K = 2L = 21 + D — 3. Still, in 
one dimension, this was just the core of feasibility of the Bender's and Boettcher's surprising 
construction. This immediately implies that the Bender's and Boettcher's one-dimensional potential 
is just a very special case of the general class of quasi-exactly solvable quartic models. Their 
description has been offered here. 

It is probably worth re-emphasizing that the vanishing of the coefficient Ak = does occur 
for the s— waves in three dimensions as well as for p— waves in one dimension. This is our re- 
interpretation of the known K = results. In the same sense, the present new K = 1 case covers 
the s— waves in four dimensions and p— waves in two dimensions. Similarly, we have K = 2 for 
s— waves in five dimensions and for p— waves in three dimensions, and we encounter the triple 
possibility of (£,D) = (0,6), (1,4) and (2,2) at K = 3. Etc. The first problem formulated in 



11 



our introductory section is satisfactorily settled: In principle, the multiplet solvability occurs in all 
dimensions. 

Our second initial motivation concerned the possible tractability of strong singularities of the 
centrifugal and Coulombic type. We have preserved their certain constrained variability within the 
Buslaev's and Grecchi's VT regularization scheme. In this way we achieved a satisfactory balance 
in the picture given at the end of section 1. 

Further lessons from our constructdion are not quite clear. In the future, we intend to pay more 
attention to the underlying Fourier-mediated p <-> x symmetries, trying to move beyond their known 
harmonic-oscillator, BG-transformation [2] or quasi-harmonic-oscillator [19] exemplifications. In 
such a context, the role of the centrifugal-like singularities does not seem to have said its last word 
yet. 

Acknowledgements 

My thanks belong to Rajkumar Roychoudhury (ISI Calcutta) and to Francesco Cannata (INFN 
Bologna) who insisted that the Bender-Boettcher model deserves a deeper study. Partial support 
by the GA AS CR grant Nr. A 1048004 is acknowledged. 



12 



References 

[1] Bender C M and Bocttchcr S 1998 J. Phys. A: Math. Gen. 31 L273 
[2] Buslaev V and Grcchi V 1993 J. Phys. A: Math. Gen. 26 5541 
[3] Bender C M and Boettchcr S 1998 Phys. Rev. Lett. 24 5243 
[4] Flcssas G P and Watt A 1981 J. Phys. A: Math. Gen. 14 L315 
[5] Sotona M and Zofka J 1974 Phys. Rev. C 10 2646 
[6] Znojil M 1999 J. Math. Chcm. 26 157 

[7] Panja M, Dutt R and Varshni Y P 1990 Phys. Rev. A 42 106 

[8] Ushveridze A G 1994 Quasi-exactly solvable models in quantum mechanics (Bristol: IOPP) 

[9] Bender C M, Boettchcr S and Meisinger P N 1999 J. Math. Phys. 40 2201 

[10] Fernandez F M, Guardiola R, Ros J and Znojil M 1999 J. Phys. A: Math. Gen. 32 3105; 
Bagchi B and Roychoudhury R 2000 J. Phys. A: Math. Gen. 33 LI; 
Znojil M 2000 J. Phys. A: Math. Gen. 33 L61 

[11] Znojil M 1999 Phys. Lett. A 259 220 

[12] Znojil M 1999 J. Phys. A: Math. Gen. 32 4563; 

Znojil M 1999 arXiv quant-ph/9912027 and arXiv quant-ph/99 12079 

[13] Alvarez G 1995 J. Phys. A: Math. Gen. 27 4589; 

Cannata F, Junker G and Trost J 1998 Phys. Lett. A 246 219; 
Delabaere F and Pham F 1998 Phys. Lett. A 250 25 

[14] Znojil 1994 J. Phys. A: Math. Gen. 27 4945 

[15] Flcssas G P 1981 Phys. Lett. A 83 121 

[16] Hautot A 1972 Phys. Lett. A 38 305 

[17] Singh V, Biswas S N and Data K 1978 Phys. Rev. D 18 1901 
[18] Turbiner A V 1988 Comm. Math. Phys. 118 467 
[19] Znojil M 1981 Phys. Rev. D 24 903. 



13 



Figure 1. Integration 
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